Adhesion clusters under shared hnear loading: a stochastic analysis 
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We study the cooperative rupture of multiple adhesion bonds under shared linear loading. Simula- 
tions of the appropriate Master equation are compared with numerical integration of a rate equation 
for the mean number of bonds and its scaling analysis. In general, force-accelerated rupture is rather 
, abrupt. For small clusters and slow loading, large fluctuations occur regarding the timepoint of final 

■ rupture, but not the typical shape of the rupture trajectory. For vanishing rebinding, our numerical 
' results confirm three scaling regimes predicted before for cluster lifetime as a function of loading 

^ , rate. For finite rebinding, the intermediate loading regime becomes irrelevant, and a sequence of 

. two new scaling laws can be identified in the slow loading regime. 

^ ■ I. INTRODUCTION 

(N : 

' Cell adhesion is based on a large variety of different adhesion molecules, each of which is optimised for its specific 
biological function 1]. Most adhesion molecules have evolved to operate under force. For example, in cell-matrix 
Q ' adhesion, receptors from the integrin family usually function under conditions of cellular contractility j^, while in 
c/2 ' leukocyte adhesion to blood vessel walls, receptors from the selectin family operate under shear flow j^]- During 
recent years, single molecule force spectroscopy has revealed formerly hidden properties of many different adhesion 
molecules, which in the future might be linked explicitly to their biological function 4]. Rupture of molecular bonds 
under force is a stochastic process which can be modeled with Kramers theory as thermally activated escape over a 
sequence of transition state barriers The most convenient loading protocol is a linear ramp of force, both 

, experimentally and theoretically. For single molecules, the most frequent rupture force as a function of the logarithm 
Q of loading rate has been predicted to be a sequence of linear parts, each of which corresponding to one transition state 
O : barrier along the rupture path 0. This prediction has been confirmed experimentally for many different adhesion 
~'i systems 0,0, including as /3i -integrin and L-selectin mediated bonds [lol| . 

I ' Although single molecule force spectroscopy has strongly changed our understanding of specific adhesion, cell 
^ . adhesion is usually not based on single molecules, but on clusters of varying size. Therefore, future understanding 
CN| ' of cell adhesion also has to include the cooperative behaviour of adhesion molecules under force. In single molecule 
, experiments, ruptured bonds usually cannot rebind due to elastic recoil of the transducer. In contrast, ruptured 

■ bonds in adhesion clusters can rebind as long as other bonds are still closed, thus holding ligands and receptors in 
' close proximity. For adhesion clusters under constant loading, it is well known that despite rebinding, stability is lost 
' beyond a critical force dj . For adhesion clusters under linear loading, force grows without bounds and the cluster will 
I always rupture. Recently, the most frequent rupture force has been measured as a function of loading rate for clusters 
. of ay/33-integrins and RGD-lipopeptides loaded through a soft transducer in a homogeneous way Theoretically, 
' it has been shown before that different scaling regimes exist for cluster lifetime as a function of cluster size, loading 

^ , rate and rebinding rate For the case of a stiff transducer, force on single bonds is independent of the number 
I ' of closed bonds and a mean field approximation can be applied to make further theoretical progress U3|. However, 
' O _ for the case of a soft transducer, force is shared between closed bonds, leading to real cooperativity: if one of the 
closed bonds ruptures, force is redistributed over the remaining closed ones. Here we present for the first time a 
full treatment of this case. We start with a one-step Master equation with Kramers-like rates, which is solved by 
Monte Carlo methods. These results are then compared to numerical integration of a rate equation for the mean 
number of bonds. We show that considerable differences exist between the stochastic and deterministic treatments 
for small clusters or slow loading. For the case of vanishing rebinding, our results confirm the three scaling regimes 
5^ ' for cluster lifetime as a function of loading rate, which have been predicted before on the basis of a scaling analysis of 
5^ I the rate equation for the mean number of bonds . For the case of finite rebinding, the intermediate scaling regime 
becomes irrelevant. For slow loading, we identify a sequence of two new scaling laws, which result from stochastic 
decay towards an absorbing boundary and finite rupture strength at constant loading, respectively. 



II. MODEL 



We consider a cluster with Nt parallel bonds. At any time t, i bonds are closed and Nt — i bonds are open 
(0 < i < Nt). The i closed bonds are assumed to share force F equally, that is each closed bond is subject to the 
force F/i. In the following, we will consider linear loading, that is F = rt where r is loading rate. Single closed bonds 
are assumed to rupture with the dissociation rate k = k^e^^^^'' , which corresponds to the case of one sharp transition 
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FIG. 1: (a) Mean number of closed bonds A'^ as a function of time r for the case of vanishing rebinding, 7 = 0, for ^/Nq = 0.1 
and Nt = A^o ~ 2, 10, lO'^ and 10''. Dotted lines: Initial exponential decay. Dashed lines: first moment of the Monte Carlo 
simulations. Solid lines: numerical integration of deterministic rate equation (not for Nt = 2). (b) Individual trajectories from 
Monte Carlo simulations in comparison with the mean (not for Nt = 2). 



state barrier along the rupture path Here, Ft is the internal force scale of the bond set by the barrier. Single 

open bonds are assumed to rebind with the force independent association rate kon- We now introduce dimensionless 
variables: dimensionless time r = fegt, dimensionless loading rate = r/k^Fb and dimensionless rebinding rate 
7 = kon/ko- The stochastic dynamics of our model is described by a one-step Master equation 

-r- = n+iPi+i + gi-iPi-i - [ri + gi\pi (1) 

CtT 

where Pi{T) is the probability that i closed bonds are present at time r. The reverse and forward rates between the 
different states i follow from the single molecule rates as 

n = ie''^/* and = 7(A^f - i) . (2) 

For constant force, this Master equation has been studied before [T3 . IT5| . Since adhesion clusters (like single molecules) 
usually cannot rebind from the completely dissociated state due to elastic recoil of the transducer, we implement an 
absorbing boundary at i = by setting go = 0. Since force increases in time without bounds, the cluster will always 
dissociate in the long run, that is Piir) — > for t oo, both for absorbing and reflecting boundaries. Cluster lifetime 
T is the mean time to reach the absorbing state i = 0. By defining cluster dissociation rate D — dpo/dr — ripi, cluster 
lifetime follows as T = Jp°° dr tD. Since the reverse rates ri are non- linear in i and time-dependent, an analytical 
solution for the pi as a function of the three model parameters A^f, /i and 7 seems to be impossible. Therefore we 
solve the Master equation numerically using the Gillespie algorithm for efficient Monte Carlo simulations, typically 
averaging over 10^ simulation trajectories for each set of parameters |l6l |. 

A quantity of large interest is the mean number of closed bonds, N = {i) = X^iJi^Pi- ^ continuum approach, 
one expects that this quantity satisfies the ordinary differential equation 

^ ^ -iVe^"/^ + 7(iVt - N) . (3) 
dr 

Cluster lifetime T can be defined by N(T) = 1. Several different scaling regimes for T as a function of Nt, fj, and 7 
have been predicted on the basis of Eq. ^ ■ Below these scaling predictions will be compared to both numerical 
integration of the deterministic equation and to our stochastic results. 



III. DECAY WITHOUT REBINDING 



We first consider the case of vanishing rebinding, 7 = 0. In this case the total number of bonds Nt does not appear 
in the model equations and the initial condition A^(0) = No is the only relevant parameter concerning the number of 
bonds. The scaling analysis of Eq. (PJ suggests that decay can be divided into two parts 0|. Initial decay is not yet 
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FIG. 2: Solid lines: deterministic results for (a) cluster lifetime T and (b) rupture force F = /iT for the case of vanishing 
rebinding, 7 = 0, as a function of ^/Na for Nt = No = 10, 10^, 10"^ and 10*. Dashed lines: curves for all three scaling regimes. 

affected by loading and thus is exponential with N[t) = Nqc^'^ . The second part of the decay is super-exponential 
and can be shown to be much shorter than the first one. Therefore the crossover time, which is defined by an implicit 
function, determines cluster lifetime T. In the regime of slow loading, /i < 1, exponential decay persists until N(t) = 1 
and T = IuNq. In the regime of intermediate loading, 1 < < A^o, the crossover occurs before N{t) = 1 is reached, 
and lifetime is reduced to T ^ ln(A'o//i). In the regime of fast loading, fi > Nq, lifetime scales even stronger with 
loading rate, T - {No/fi) \n{fi/No). 

In Fig. ^ we plot N{t) as obtained from simulations of the Master equation (dashed lines) and from numerical 
integration of the deterministic equation (solid lines) for Nq — 2, 10, 10^ and lO'^. The dotted lines are the exponential 
decays N{t) — NQe~'^ for vanishing loading. In the presence of loading, the later part of the decay process clearly 
is super-exponential. The first moment of the stochastic process decays less abrupt than the deterministic result, 
although for increasing cluster size, the difference between stochastic and deterministic results becomes smaller. In 
Fig. we show representative trajectories from Monte Carlo simulations. They demonstrate that the final stage of 
the rupture process is rather abrupt. In fact abrupt decay is typical for shared loading and is found also for shared 
constant loading: a decreasing number of closed bonds increases force on the remaining bonds, thus further increasing 
their dissociation rates jl5l |. As Fig. ^ shows, fluctuations tend to change the timepoint of rupture, rather than 
the typical shape of the decay curve. For increasing cluster size, fluctuations become smaller and rupture events are 
concentrated around the rupture of the deterministic cluster. An analysis of the variance of the number of closed bonds 
i shows that for slow loading, it is close to the exact result for vanishing loading, (i^) — (i)^ = Noe^'^{l — e^"^) [l^ . 
It vanishes for r = due to the initial condition, then quickly rises to a maximum and finally decays exponentially. 
As loading rate fi increases, a large additional peak appears shortly before final rupture (not shown). 

Simulations allow to measure cluster dissociation rate D{t) and cluster lifetime T for all parameter values. For 
/i < 1, the simulation results are close to the known analytical results for /x = 0, D{t) — Nqc^'^^I — e~'^)'^°~^ and 
T = X;,=°i 1/i « In A^o + (l/2iVo) -h 0.577 (1^,19). For large A^o, the deterministic scaling T^IuNq results. For fi>l, 
the functions D{t) become narrowly peaked around the mean value T. As suggested by the scaling analysis, we find 
that now T depends only on the value of /i/iVo. In Fig. we plot deterministic results for T as a function of /i/A^o 
and for different values of Aq. The stochastic results are very similar, except for the differences in the initial plateau 
values. Initially, the different curves plateau at the values InAo for /i < 1. For 1 < fi < Nq and sufficiently large Aq, 
they collapse onto a universal curve, which can be approximated by 0.841n(0.35A'o//i). For /i > Nq, they collapse 
onto another universal curve, {Nq / ^) l'!i{fi / Nq) . In Fig.|2j3, we plot the logarithm of the deterministic rupture force, 
F — jiT, as a function of h/Nq. For large Nq, one clearly sees the sequence of the three different scaling regimes. For 
decreasing A'oi the intermediate scaling curve becomes an increasingly bad fit. 

IV. EFFECT OF REBINDING 

In the stochastic framework, cluster lifetime can be identified with the finite mean first passage time of reaching 
the absorbing boundary at i = 0. For ^ = Q and A'o = Nt, an exact result can be obtained with the help of Laplace 
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FIG. 3: (a) Mean number of closed bonds A'' as a function of time r for rebinding rate 'y — 1 and loading rate /i/No — 0.01. 
Cluster sizes Nt = 2, 10, 10^ and 10'^, initial condition A'o = Neq — ^Nt/{1 + 7). Dotted lines: Initial number of closed bonds. 
Dashed lines: first moment of the Monte Carlo simulations. Solid lines: numerical integration of deterministic equation (not 
for Nt = 2). Dashed-dotted lines: effect of reflecting boundary for Nt — 2 and 10. (b) Individual trajectories from Monte Carlo 
simulations in comparison with the mean (not for Nt =2). 
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where Hfft = X]iJi(V*) is the iVt th harmonic number. In the deterministic framework of Eq. ^ and /it = 0, an 
adhesion cluster with a total of Nt molecular bonds will equilibrate from any initial number of closed bonds A^o to a 
stable steady state with N^q — jNt/{l + 7) closed bonds. For convenience, in the following we will use A'o = N^q. 
Then similar results follow for Tstoch as given in Eq. A stability analysis of the deterministic equation Eq. (jSj 
for loading with a constant force f = F/Fh shows that the steady state cluster size decreases until stability is lost 
beyond a critical force fc — A^tplog(7/e) [lllial- Here the product logarithm plog(a) is defined as the solution x of 
xe^ = a. For 7 < 1, the critical force can be approximated as fc ~ Nt"f/e: it vanishes with 7 since without rebinding 
the cluster decays by itself. For 7 > 1, it can be approximated as fc ~ 0.5Nt In 7, that is the critical force now is only 
a weak function of rebinding. For slow loading, /i < 1, the adhesion cluster will follow the quasi-steady state until the 
critical force fc is reached at the time = /c//^- The remaining time to rupture is smaller and thus the lifetime of 
the adhesion cluster is close to 
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It diverges with the inverse of loading rate in the limit of vanishing /it, as it is required by the existence of a stable 
steady state and a finite rupture force F = fiT = fc. For intermediate loading, 1 < < TVq, a power-law behaviour 
T ^ {Nq/h)^^^ has been erroneously predicted in ref. as reported in ref. 13]. For fast loading, fi > Nq, rebinding 
can be neglected and T ~ {No/ 11) ln(/z/A^o) as in the previous section. 

In Fig. 12^ N{t) is plotted for 7 = 1 as obtained from Monte Carlo simulations (dashed lines) and from numerical 
integration of the deterministic equation (solid lines). The different initial conditions Nq for Nt = 2,10,10'^ and 
10'^ are represented by the dotted lines. In Fig. ^jp individual trajectories from the simulations are compared to the 
stochastic averages from Fig. For the small clusters, Nt = 2 and 10, loading rate is so small that fc/fJ- > Tstoch- 
Then T « Tstoch and the clusters decay by themselves due to stochastic fluctuations to the absorbing boundary (ultra- 
slow regime). The dash-dotted lines in Fig. |3l show the effect of a reflecting boundary, which is rather dramatic for 
these small cluster sizes. For the large clusters, Nt = 10^ and lO'^, fluctuations are less probable until the force is 
close to fc. Therefore the individual clusters fluctuate around the quasi-steady state and dissociate only close to the 
deterministic cluster lifetime Tdet. Due to the large force on a single bond at fc, the boundary has little influence 
here. A detailed analysis of the variance of i confirms this description (not shown): for the smallest cluster, when 
fluctuations dominate during the whole time evolution, the variance shows a broad peak. For the larger clusters, it 
develops a narrow peak around the mean rupture time. 
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FIG. 4: (a) Mean cluster lifetime T and (b) mean rupture force F — jiT for the case 7 = 1 as a function of fi/No for 
Nt = 2, 10, 10^ and 10^ In (a), the curves for the two larger clusters are nearly identical. 

In Fig. 21 we show mean cluster lifetime T and mean rupture force F — /iT for 7 = 1 as a function of ^/Nq for the 
cases Nf = 2, 10, 10^ and 10"^. For the small clusters, T starts at the value of Tgtoch and ends in the scaling regime 
for fast loading, where the curves are practically identical for all different parameter values at a given value for (i/Nq. 
The curves for the large clusters are nearly identical. They start at the values of Tdet for small loading rates and 
end in the same fast loading regime. An intermediate loading regime seems to exist only as a transient between the 
regimes of slow and fast loading. In particular, it does not fit well to an inverse square root dependence, as shown in 
Fig. Hi. 



CONCLUSIONS 



In this paper, we have presented for the first time a full analysis of the cooperative decay of a cluster of adhesion 
bonds under linearly rising force. Significant differences between stochastic and deterministic treatments are found 
for small clusters or slow loading, when stochastic fluctuations are relevant. However, they do not affect so much the 
typical shape of the rupture trajectory, but rather the timepoint at which rupture occurs. For the case of vanishing 
rebinding, 7 = 0, our full treatment nicely confirms the scaling analysis of the deterministic equation for cluster 
lifetime T as a function of fi and A^o |l2l |. However, in contrast to the scaling analysis, the full treatment presented 
here allows for detailed comparision with experiments, e.g. in regard to typical unbinding trajectories or binding 
strength over a range of loading rates spanning different scaling regimes. For the case with finite rebinding, 7 > 0, 
we identify a sequence of two new scaling laws within the regime of slow loading, fj, < 1. For ultra-slow loading, T is 
independent of /i and is determined by stochastic fluctuations towards the absorbing boundary. For larger fi (but still 
with < 1), T starts to scale inversely with /x, due to the flnite rupture strength at constant loading. In contrast to 
the case of vanishing rebinding, a scaling regime of intermediate loading, 1 < /x < Nq, could not be identified. 

Our results can be applied for example to rolling adhesion of leukocytes, when multiple L-selectin bonds are 
dynamically loaded in shear fiow |20| |. Dynamic force spectroscopy has only recently been applied to clusters of 
adhesion bonds 113. RGD-lipopeptides on a vesicle have been presented to a^/Ja-integrins on a cell. The effect of 
thermal membrane fluctuations can be disregarded on both sides, because the vesicle is under large tension and the 
intcgrins are rigidly connected to the cytoskeleton. Appreciable loading occurs only over a ring region along the rim 
of the contact disc, for which no inhomogeneities have been observed. If one neglects the subsequent peeling of the 
inner region, which presumably is much faster, our model can be applied. The parameter values can be estimated 
to be Nt « 100, Ff, « 40 pN, ko « 0.01 Hz and 7 w 1. Loading rates have been varied from r = 20 — 4 x 10'^ 
pN/s, that is fJ./Nt — 0.5 — 100. Therefore this experiment should correspond to the intermediate and fast loading 
regimes. We expect that future improvements in experimentation will make it possible to probe also the slow loading 
regime, where rebinding and stochastic effects become relevant. In order to achieve a more complete understanding 
of the role of force in cell adhesion, future modeling should also address the detailed nature of the force transducer, 
non-homogeneous loading and more realistic scenarios for the rebinding process. 
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